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Abstract 

We study the long-distance dominated Kl — > 7r"'"7r~7 decay at 0{p^) in Chiral 
Perturbation Theory. The complete calculation of the 0{p^) loop magnetic am- 
plitude is carried out. At this chiral order the model dependent part of the vector 
meson exchange contribution to the magnetic amplitude is evaluated in the two mod- 
els : FM (Factorization Model) and FMV (Factorization Model in the Vector Cou- 
plings). We predict, in an almost model independent way, a slope of Kl 7r+7r~7 
in the range c ~ —1.6,-1.8, consistently with the experimental value. We find that 
the experimental result for the width of Kl tt~^tt~'j is compatible with a bigger 
breaking of the nonet symmetry in the weak vertices than previously stated. Thus 
we conclude that our analysis does not exclude an opposite sign to the one given by 
the dominance of the pion pole in the poorly known Kl — > 77 amplitude. A com- 
plete analysis of the 0{p^) weak Vector-Pseudoscalar-Pseudoscalar (VPP) vertex 
is also performed. 
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1 Introduction 



Radiative non-leptonic kaon decays provide relevant tests for the ability of Chiral 
Perturbation Theory (xPT ) [|l], Q to explain weak low-energy processes. xPT is a natural 
framework that embodies together an effective theory, satisfying the basic chiral symmetry 
of QCD, and a perturbative expansion in masses and external momenta that becomes the 
practical tool to work with. Its success in the study of radiative non-leptonic kaon decays 
has been remarkable (see Refs. 0, |] and references therein). 

Moreover the chiral anomaly present in the Standard Model manifests mainly in the 
low-energy interactions of the pseudoscalar mesons, and consequently xPT is the appro- 
priate framework to study those effects. As has already been shown 0, ^ , radiative kaon 
decays are sensitive to the chiral anomaly in the non-leptonic sector. The study of these 
decays driven by the chiral anomaly {K mr'-f, K mi'-f'j, etc.) has thoroughly been 
carried out previously in the Refs. [@, ^ |3, where their relevant features have been 
discussed. 

The main uncertainty still present in the computation of those channels is the contri- 
bution from local terms that, in principle, appear at any order in the chiral expansion. At 
0{p^) vector meson exchange (when present) has been shown to be effective in predicting 
the relevant couplings in the strong sector However our poor phenomenological 

knowledge of the weak processes involving meson resonances translates into our ignorance 
about the couplings in weak counterterms at 0{p'^) and beyond, thus models are required. 
In particular the Factorization Model (FM) has been widely employed in this task 



|12| , [T3|. In Ref. |^ we proposed an implementation of the FM in the Vector Couplings 
(FMV) as a more efficient way of including the contributions of vector mesons into the 
weak couplings at 0{p^) in the processes K n'yy and Kl — > ■y£~^£~. We reached a 
good phenomenological description of both channels and showed that with the new vec- 
tor contributions included no enhancement due to AI = 1/2 transitions is required, in 
opposition to previous statements in the FM. 

The basic statement of the FMV is to use the idea of factorization in order to construct 
the weak vertices involving vector mesons instead of implementing factorization in the 
strong lagrangian generated by vector exchange once the vector degrees of freedom have 
already been integrated out. The apparently tiny difference between both realizations 
has proven to be crucial in disentangling the role of the FM in uncovering new chiral 
structures that were absent in the standard approach and consequently giving a deeper 
understanding of the underlying physics. 

In the present article we present the application of these ideas to the process 
7r+7r^7 . This decay brings new features to our previous work. While in ^ 7177 and 
Kl —>■ 'ji^i" only the weak 0{p^) Vector-Pseudoscalar-Photon {VP'j) vertex is involved, 
now in Kl —>■ 7r+7r^7 we also need to consider the weak 0{p^) Vector-Pseudoscalar- 
Pseudoscalar {VPP) vertex. The interest in studying this particular process has a twofold 
purpose: a) the relevance of the process by itself given the fact that phenomenologically 



constitutes an appropriate target of xPT , and b) the dependence of this channel in the 
breaking of nonet symmetry in the weak sector. Let us comment these two aspects in 
turn. 

As studied time ago |]l5l the process ti^h^^ is driven by long-distance contri- 

butions. The Direct Emission amplitude of the decay is dominated by a magnetic type 
amplitude which slope in the photon energy in the kaon rest frame (to be defined properly 
in Section 4) has been measured experimentally. Even when the error is still a rough 30 % 
the slope seems rather big, a feature to be explained in xPT as a contribution starting 
at 0{p^) . Hence this fact together with our poor knowledge of the constant amplitude 
deserves a careful treatment. At 0{p^) chiral loops are present too, thus we have also 
performed the full loop evaluation at this chiral order. It turns out to be relevant for the 
stability of our slope prediction. 

In Ref. fl^ we pointed out that the experimental slope of Kl 77* and the theoret- 
ical vector meson dominance prediction cannot be accommodated unless the Kl 77 
amplitude departs drastically (change sign) from the one predicted by the nonet symme- 
try. We would like to suggest that at the origin of this problem could be a large breaking 
of the nonet symmetry in the weak vertex. A way to explore this possibility is to study 
this hypothesis in another process sensitive to this breaking: Kl 7r+7r~7 receives an 
0{p^) contribution due to anomalous reducible amplitudes f^, |^ where the nonet break- 
ing plays a major role. Therefore we would like to investigate also the possibility of a 
satisfactory simultaneous description of the processes Kl — > 77 and Kl — > 7r+7r^7 . 

In Section 2 we will remind briefly the basic features of xPT and its application to 
weak processes. Then in Section 3 we will specify the general characteristics and notation 
for Kl 7r+7r^7 and we will offer a short overview of the process in xPT . The 0{p^) 
contributions to the magnetic amplitude will be collected and explained in Section 4. 
The analysis of the combined observables will be carried out in Section 5 while we will 
emphasize our conclusions in Section 6. Three brief appendices complement the main 
text. 



2 Non-leptonic weak interactions at low energies 

We will review in this Section the procedures of xPT and their implementation in the 
study of non-leptonic weak processes. We will collect also the tools we will need in the 
development of our study. 

2.1 Chiral Perturbation Theory 

XPT m, ^ is an effective quantum field theory for the study of low energy strong 
interacting processes that relies in the exact global chiral symmetry of massless QCD. 
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The basic assumption is that, at low energies {E < 1 GeV), the chiral symmetry group 
G = SU{3)L^SU{3)ii is spontaneously broken to the vector subgroup SU{3)v, generating 
eight Goldstone bosons to be identified with the lightest octet of pseudoscalar mesons {(pi). 
These are conveniently parameterized by a SU (3) matrix field 



U{ip) = exp - Y.^j'Pj ' (1) 




that transforms linearly under the chiral group G. In Eq. (|l]) Aj are the SU{3) Gell-Mann 
matrices Q and F ~ ~ 93 MeV is the decay constant of the pion. 

The extension of the global chiral symmetry to a local one and the inclusion of external 
fields are convenient tools in order to work out systematically the Green functions of 
interest and the construction of chiral operators in the presence of symmetry breaking 
terms. A covariant derivative on the U field is then defined as 

D^U = d^U - ir^U + iUi^ , (2) 

where = — and = f/^ + a^i, are the left and right external fields, respectively, 
in terms of the external vector and axial-vector fields. If only the electromagnetic field 
is considered then = = —eQA^ where Q = diag{2/3, —1/3, —1/3) is the electric 
charge matrix of the u, d and s quarks 0. The explicit breaking of the chiral symmetry 
due to the masses of the octet of pseudoscalars is included through the external scalar 
field s = M. + ... \\i this way an effective lagrangian can be constructed as an expansion 
in the external momenta (derivatives of the Goldstone fields) and masses JI], 0, |I6[. 
The leading OijP') strong lagrangian is 

A = -^{uX + X+) , (3) 
where [A) = Tr{A) in the fiavour space and 
= iv) U It 

X+ = Xu^ + ux^ u , X = '2Bo{s + ip) = 2BoM + ... , (4) 

M = diag{mu , , m^) , i?o = - ^ (0| uu |0) . 

r 

The even-intrinsic parity lagrangian at was developed in Ref. and introduces 12 

new coupling constants. For our work the only piece we will need is 

A = -iLg{f1''Uf,Uy) + ... , (5) 

where 

/r = uF^^u^ ± u^F^^u , (6) 



1 Normalized to Tr(\i\j) = 2Sij. 

^This corresponds to Df^cj)^ — (9^ ± ze^^)(^^. 
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and F^j^ are the strength field tensors associated to the external and £^ fields, re- 
spectively. From the experimental value of the pion charge radius one obtains Lg(mp) = 
(6.9 ±0.7) X 10-3. 

The 0{p'^) odd-intrinsic parity lagrangian arises as a solution to the Ward condition 
imposed by the chiral anomaly [|17|. The chiral anomalous functional Za„[f/, r] as given 
by the Wess-Zumino-Witten action (WZW) is 



z407r JM^ 



(7) 



4871^ 



W{U,i 

-(L^i?), (8) 

with Nc = 3, E^ = Wd^U, E^ = Ud^W and (L ^ R) stands for the interchange 
U <-> W, ^ r^, E^ ^ E^. We notice that the WZW action does not include any 
unknown coupling. 

The inclusion of other quantum fields than the pseudoscalar Goldstone bosons in the 



chiral lagrangian was considered in Ref. |1T8|. We are interested in the introduction of 



vector mesons coupled to the U{ip) and to the external fields. Let us then introduce the 
nonet of vector fields 



1 A. ... 1 



that transforms homogeneously under the chiral group G. Here ideal mixing, i.e. V"^ = 
{u^ + \/2(f)^) / y/S, is assumed. 

The most general strong and electromagnetic lagrangian, at leading 0{p^) and assum- 
ing nonet symmetry, with the vector field linearly coupled to the Goldstone bosons, is 



given by the following terms |19 



+ ieve^,p^{V''u''u''u'') , (10) 



where only the relevant pieces for our work have been written. In Eq. (10) V^^, = V^VJy 
VjyV^ and is the covariant derivative defined in Ref. |^ as 
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(11) 

T/, = - { u\d^ - ir^)u + - if^)^"^ } . 

The couplings in Eq. ([T0|) can be determined phenomenologically from the experiment pO[| . 
The experimental width T{uj tt^'-^) gives \hv\ — 0.037 0]. Also from r(p'^ — e+e^) 
one gets — 0.20 and from r(p° ^ tt+tt^) l^fyj ~ 0.09 0. Moreover the positive 
slope of the tt^ —>■ 77* form factor determined experimentally imposes hyfv > 0. The 
resonance saturation of the 0{p^) couplings gives |^ Lg = fvgvf^ in Eq- (H) and the 
positivity, determined phenomenologically, of U^ijUp) implies fvgv > 0. In the Hidden 



Gauge Model the relations fy = 2gv and Oy = 2hv arise. The first one is in excellent 
agreement with the experimental determination while the second is in good agreement 
with the ENJL model results We will assume them throughout the paper. 



The introduction of the axial-vector nonet defined analogously to the vector in Eq. 
is similar to the latter. The interaction terms we will need are 



V/-' ) + hA e,..p ( {u\ r}), (12) 



with A^i, defined as in the vector case above. The couplings /a and Ha could be determined 
from radiative decays. Due to the poor phenomenology available on axial-vector decays 
we rely on theoretical predictions |2^ that give ~ 0.09 and — 0.014. 



The incorporation of spin-1 mesons in chiral lagrangians is not unique and several 
realizations of the vector field can be employed [^. In Ref. [0 was proven that at 



in xPT and once high energy QCD constraints are taken into account, the usual 
realizations (antisymmetric tensor, vector field, Yang-Mills and Hidden formulations) are 
equivalent. Although the antisymmetric tensor formulation of spin-1 mesons was proven 
to have a better high-energy behaviour than the vector field realization at 0{p'^) , this fact 
is not necessarily the case in general. In fact for the odd-intrinsic parity operator relevant 
in the V — > P7 decay, the antisymmetric tensor formulation would give contributions 
starting at O(p^) while QCD requires an explicit 0{p^) term as given by the vector 
realization in the hy term in Eq. ( p!OD pT| Moreover we have pointed out that the 
conventional vector formulation seems to give a complete and consistent treatment of the 
spin-1 resonance generated weak couplings at 0{p^) p6[ and 0{p^) |l^ . 
For a further extensive and thorough exposition on xPT see Ref. [p7[]. 



^For a more detailed discussion of the equivalence of vector resonance models in the odd-intrinsic 
parity violating sector see Ref. ||2^. 
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2.2 Non-leptonic weak interactions in 

At low energies {E <^ M^y) the AS* = 1 effective hamiltonian is obtained from the 
Lagrangian of the Standard Model by using the asymptotic freedom property of QCD in 
order to integrate out the fields with heavy masses down to scales /i < rric. It reads 

n^N^' = K.V;; E + h.c. . (13) 

Here Gi;' is the Fermi constant, Vij are elements of the CKM matrix and Cj(/i) are the 
Wilson coefficients of the four-quark operators Qi, = 1, •••6. ?^1vl induces AJ = 1/2 
and A/ = 3/2 transitions. As the first components are phenomenologically dominant 
in the non-leptonic kaon decays and since data in the structure dependent amplitude of 
Kl —>■ 7r+7r~7 are accurate only until 20% level, the AJ = 3/2 transitions can be neglected 
for the present purposes. Therefore we only will consider the octet component of 

At 0{p'^) in xPT we can construct only one relevant octet effective operator using the 
left-handed currents associated to the chiral transformations, 

4^""=' = 4^V^,V:,gs{XeLlLl) 

= ^F'V^,V:,gs{AuX) , (14) 

where 

Ll = J^ = -^^U^D.U = -—u^u.u , (15) 

is the left-handed current associated to the 5*2 action of the 0{p'^) strong lagrangian in 
Eq. (H), A = uXqu^ and is an effective coupling. 

From the experimental width of — > vrvr one gets at C(p2) 

l^sk^^vr ^5.1 , G^ = ^KdK; \gs\K^^. ^ 9.2 x 10"^ GeV-^ . (16) 

In our analysis of 7r"^7r~7 the singlet pseudoscalar r/o contributes at 0{p^) . If nonet 

symmetry is broken weak interactions of the singlet at 0{p'^) have to be parameterized 
through a new coupling : p. Then 

^^'^'^^ = \G,F'{p-l){/\u,){un . (17) 
should be added to £3^'^'"^ in Eq. ([l 



At 0{p ) the weak chiral lagrangian has been studied in Refs. W\ B^, RD[ giving 37 



chiral operators only in the octet part. For the study of anomalous radiative decays only 
four of them are relevant. 



+ Nso{Aun{ffun + iV3i(Aw^)(/^V) 



+ ... 
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where A^28) ••• 5-^31 cire new coupling constants to be determined phenomeno logically or 
predicted by models. 

Motivated by l/N^. arguments the concept of factorization has been developed in the 
context of xPT . In the Nc ^ 00 limit penguin operators are suppressed. Then if the 
effect of gluon exchange at leading order is considered and we assume octet dominance 
through current-current operators in Eq. (|l^) we have ||2^ 



= - ^ v-^v:s c-{ii) Q- + h.c. , (19) 

where 

g_ = 4 ( sli'^ul ) ( ULlf.dL ) - 4 ( SLl'^di ) ( mlTm^l ) , (20) 

with qiL'~i^q2L = |?i 7^(1 "75)920 and a a colour index. The Q- operator transforms under 
SU{3)l SU{3)fi as the (8^,1^?) representation. The result for the Wilson coefficient 
C-{mp) at leading 0{as) and with A;^ = 325 MeV is 

C_(m,) ^ 2.2 gs-l.l. (21) 

to be compared with Eq. ([T6|). In a chiral gauge theory the quark bilinears in the Q- 
operator in Eq. (|^) are given by the associated left-handed current 

§ = L; + Lj + i^+-. (22) 

(the first term was already given in Eq. ([T5|)) where ^[f/, £, r, is the low-energy 
strong effective action of QCD in terms of the Goldstone bosons realization U and the 
external fields r, The assumption of Factorization amounts to write the four-quark 
operators in the factorized current x current form as 

= 4A:^G8(A^^) + h.c. , (23) 

where A = ^(Ae — iXj) and /cj? is an overall factor not given by the model. In general 
kp — 0{1) and naive factorization would imply fci? ~ 1. However the Wilson coefficient of 
the Q- operator would give kp — 0.2 — 0.3 as follows from Eq. (^Tf ) and thus a dynamical 
understanding of factorization is achieved only for a value close to this one. Indeed we 
will find for our factorization results this value. Therefore we imply that the A/ = 1/2 
enhancement in K ^ tttt is specific of this channel but it is not a general feature of the 
FM in all non-leptonic kaon decays. 

It has become customary in the literature p, |^ to introduce the a, couplings defined 
through 

^V28 - ^ , ^V29 - —2 , 



A'ao = , ^n = ^. (24) 
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In Refs. |3T| the FM was used to evaluate the factorizable contribution of the chiral 
anomaly (WZW action in Eq. (|^)) to these couplings, giving a"" = rjan, i = 1,2,3,4, 
with 7]an ~ and positive the unknown FM factor. Non-factorizable contributions do 
not add any new chiral structure to these. 

This same model predicts vanishing resonance contributions to these couplings if the 
antisymmetric formulation for the spin-1 fields is used. In Ref. we have evaluated 
the spin-1 resonance contributions (vector and axial-vector) to the Ni couplings using a 
novel framework in which the vector formulation of the resonance fields is implemented. 
We quote here the relevant results for our case using factorization. These are 

= V2fve 



vVv , 

N29^^ = -^ifvhvVv - Ia hA Va ) 



= 2V2 fvh 



vVv 



= 2V2 f A hAVA, (25) 



where fv, hy, Oy, /a and Ha have been defined in Eqs. (|T0|,|T2|), and rjy, rjA are the 
undetermined factorization parameters. We have, therefore, three unknown couplings : 
Van, Vv and r/A- We remind that naive factorization, however, implies = = Va = V- 
Once we determine the couplings from the phenomenology of 7r+7r~7 (see 

Section 5) for sake of definiteness we will assume the naive FM relation. In this way we 
reduce the four Ni couplings to only one free parameter. 



3 The Kl n^n 7 amplitudes in 

The general amplitude for K irn'j is given by 

A[K{p)^7rM7r,{p2Hq,e)] = e^*{q) M^{q,p,,p2) , (26) 

where e^(g) is the photon polarization and is decomposed into an electric E and a 
magnetic M amplitudes as 

= [pi ■ gp2M - P2 ■ qPl^,] + ^^3'^ Bf.upaP'ipk'' , (27) 

171}^ m]^ 

with 

Zi = — — , (z = 1, 2) , Z3 = , 23 = zi + Z2 . (28) 

"i^K rriK 

The invariant amplitudes E{zi), M{zi) are dimensionless. Adding over photon helicities 

there is no interference between both amplitudes and the double differential rate for an 

unpolarized photon is given by 

d'^T _ rriK 
dzi dz2 (Att)^ 



{\E{zi)f + \Mizi)f) [z,Z2il - 2{zi + Z2) - rl - rl) ~ rlzl - rlzl 



(29) 
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where rj = uiT^./mK- In Appendix A we recall, for completeness, some of the kinematical 
relations for tt^tt^'j . 

The total amplitude of K ^ TTTC'-f can be decomposed as a sum of inner hremsstrahlung 
(IB) and direct emission (DE) (or structure dependent) amplitudes. The electric ampli- 
tude Eisizi) arises already at C(p^) in xPT and it is completely predicted by the Low 
theorem ||3^ which relates radiative and non-radiative amplitudes in the limit of the 
photon energy going to zero. Due to the pole in the photon energy the IB amplitude gen- 
erally dominates unless the non-radiative amplitude is suppressed due to some particular 
reason. This is the case of Kl 7r+7r~7 (where K2 tt+tt^ is CP violating 0, |[) or 
7r+7r°7 (where — > tt+tt" is suppressed by the A/ = 1/2 rule). These channels 
are important in order to extract the DE amplitude that reveals the chiral structure of the 
process. In this article we are going to focus our attention in the Kl — > 7r^7r~7 channel 
(for a thorough overview of other channels see Refs. 0, §]). 

DE contributions can be decomposed in a multipole expansion |3^, In our following 



discussion we will consider only CP-conserving DE amplitudes. These start at 0{p'^) in 
xPT where E2 and Ml multipoles are generated. Due to the asymmetry under the 
interchange of pion momenta of the E2 amplitude there are no local contributions and 
then it is only generated by a finite chiral loop amplitude E^^lp. It was computed in 
Refs. I, i. The amphtude E^^l^ gives an E2 multipole very suppressed in comparison 
with the IB, typically \Ei^lp/EiB\ < 10"^. This is due to several circumstances: chiral 
loop suppression, absence of the photon energy pole and asymmetry under the interchange 
of pion momenta. Higher order terms, though in principle could be larger, still generate 
very small interference with the IB amplitude p[. Thus we will neglect here the electric 
amplitudes. 



In fact the experimental results for the branching ratio of Kl — > vr+vr 7 are 
consistent with the IB {Br{KL 7r+7T~-f ; E* > 20MeV)/B = (1.49 ± 0.08) x 10^^), and 
a dipole magnetic contribution 

Br{KL 7r+7r^7 ; E* > 20 MeV)DE = (3.19 ± 0.16) x 10"^ , (30) 

where E* is the photon energy in the kaon rest frame. 

In opposition to the electric case the Ml multipole at leading 0{p^) is generated by a 
constant local contribution only. This is given by £4^'^'"^ in Eq. ([T8|) through the diagram 
in Fig. l.a and the result is 

M(^) = -^^[32vr^(Ar,, + Ar30] 
l-n 



02 + 204). (31) 



As there is no loop contribution at this order the combination of couplings N29 + N^i 
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b) 



Figure 1: Diagrams contributing to the 0{p'^) magnetic amplitude M^^\ The diagram 
b) is vanishing at this order. The crossed box corresponds to a vertex generated by the 
^\As\-i ig^gj^g^^gjg^^ (|l8|) , thc cmptj box is generated by C^^^^~^ in Eq. (p^ and the 

circle by the WZW action Zan in Eq. (^. 



is scale independent 0. These couplings are unknown from the phenomenology and this 
fact constrains our knowledge on this contribution. At O(p^) one could also consider the 
reducible anomalous magnetic amplitude generated by the diagram in Fig. l.b. However 
at this chiral order this amplitude vanishes because the Gell-Mann-Okubo mass relation 
as happens in Kl 77 . Due to this cancellation and the possibility that the combi- 
nation 02 + 2a4 is small, the 0{p^) contributions seem theoretically important. Also the 
experimental analysis carried out in Ref. ||35[ shows a clear dependence in the photon 
energy that appears in the magnetic amplitude only at 0{p^) . It is therefore necessary 
to study the 0{p^) contributions to the magnetic amplitude in order to be able to analyse 
thoroughly this decay. This we will do in the following Sections. 



4 Magnetic amplitude of Kl tt+tt 7 at 0{p^) 

Both local terms and chiral loops contribute to a magnetic amplitude for — > 7r+7r^7 
at 0{p^) in xPT . As commented before the first dependence of the magnetic amplitude 
in the photon energy appears at this chiral order. 

*In fact all the anomalous couplings in £4^'^'"^ in Eq. ( |l8| ) are separately scale independent because 
there is no intrinsic parity violating action at 0{p^) . 
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The magnetic amplitude can be expressed as 

M = -^^m(l + c.3), (32) 

where \m\ and the slope c are fixed by the rate and the spectrum, respectively 0. From 
the experimental results in Ref. ||3^ we have 



l^lexp = 1.53 ±0.25 



Ceip — 1-7 i 0.5 , 



(33) 



where the error in \rn\f.xp comes only from the error in the slope, z^, has been defined in 
Eq. ( pSD and in the rest frame = E*/mK- 

As can be noticed the dependence on the photon energy is by no means negligible. 
However we have seen that it is vanishing at 0{p'^) in Eq. (^) and then it is a goal of the 
higher orders in xPT to accommodate such a big value. 



4.1 Local amplitudes 



There are different local contributions to Kl 7r+7r^7 starting at 0{p^) . Some of 
them are model independent while in some other cases one has to use models. Here we 
do an analysis of all of them at this chiral order. 

First of all there is a reducible anomalous amplitude given by the diagram in Fig. 2. 
This is the same that the vanishing one in Fig. l.b but now the 7r°, rj and rj' are included. 
Then it is 0{p^) (but also contains higher orders) and is given by 



(34) 



where 



1 — r„ 



[(1 



3(1 



3(1 -r; 
J- [2V2pcos^ 



Ocos0 
-(1+Osin^] 



2v^psin e\ 



F s 3 



cose - V2 



sin e 



sme + V2(^) cose 



(35) 



In Fi, rp = mp/niKi ^ parameterizes the S'f/(3) breaking defined through 









Kl) 




£|AS'|=l 


Kl) 



(36) 



^Neglecting the Electric multipoles in the total DE amplitude. 

^The value of the slope is given by E.J. Ramberg in private communication reported in Ref. Q]. The 
central value can also be followed from Ref. pal . 
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D 




7 



TT 



Figure 2: Diagram contributing to the reducible anomalous 0{p^) magnetic amplitude 
^anom- The empty box is generated by £2^'^'""^ + C^^^^~^ in Eqs. (|14|, |l^) and the white 
circle by the WZW action Zan in Eq- (0)- 

9 denotes the rj-rj' mixing angle, p has been defined in Eq. ([T7|) and Fg, Fq are the decay 
constants of rjs and rjo, respectively. At 0{p'^) Fi vanishes because the Gell-Mann-Okubo 
mass relation. We note that M^^^^ is constant and has no dependence on the photon 
energy. 

Other 0{p^) local amplitudes can be generated through resonance interchange. In 
particular we are going to focus in the vector meson contributions which are reasonably 
thought to be the most relevant ones. There are two different kinds of these contributions 



a) Vector exchange between strong/electromagnetic vertices with a weak transition in 
an external leg. These are usually called indirect amplitudes and they are model 
independent due to our good knowledge of the strong and electromagnetic processes 
involving vector mesons. 

b) The direct transitions are those where the weak vertices involving resonances are 
present and our poor knowledge of weak decays of resonances makes necessary the 
use of models in order to predict them. 

The indirect model independent vector contribution corresponds to the diagrams in Fig. 3 
that give (assuming = 0) 



■ 




27T^F^ 



rv [1 - 3^3] , 



(37) 
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Figure 3: Diagrams contributing to the indirect vector meson exchange magnetic amph- 
tude M-^l^. The empty box is generated by C^^^^~^ in Eq. ([T^) and the black circles by the 



strong/electromagnetic lagrangians Cy in Eq. (0). In b) the crossed diagram ^ n 
has to be considered too. 

rv = — ^ ~ 0.41 , 

where the numerical value of ry is for my = rrip []. The strong/electromagnetic vertices 
present in the diagrams in Fig. 3 are given by the lagrangian density Cy in Eq. (|10|). 
Several points are worth to comment: a) M-^^ shows a dependence in zs (i.e. in the 
photon energy); b) contrarily to what happened in the reducible anomalous contributions 
(Fig. l.b and Fig. 2), even considering only the 7r° and the t]s states, M-^l^ does not vanish 
and our result agrees with the one quoted in Ref. [0. 

Both contributions Mj^f}^^ and M^^\ are characterized by a weak transition in the 
external legs given by £2^'^'"^ in Eq. (|T^). We do not consider the Feynman diagrams 
where this transition happens in a pion final leg because they are suppressed by m^/m|^ 
over those where the weak vertex happens in the initial Kl leg. 

The direct resonance exchange contributions are generated by the diagrams in Fig. 4 
where a direct weak vertex involving vector mesons is present. Due to our ignorance on 



''If we include the singlet rjo and rj-i]' mixing in Fig. 3 we would have found an extra term 



~(6) GaefTiK 



(38) 



— — 3zq 

2 ^ 



to be added to Mj^^\ in Eq. (37) {Fi has been defined in Eq. (|3^)). However this term starts to contribute 



„8 



at 0{p ) and therefore is out of our scope here. 
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these vertices from the phenomenological point of view it is necessary to invoke models in 



order to predict them. Moreover it has been shown \T^, |2T| that these contributions are 
not small compared with the indirect ones and therefore have to be taken into account. 

Thus we construct the most general weak 0{p^) Vector-Pseudoscalar-Photon (VP'y) 
and 0{p^) Vector-Pseudoscalar-Pseudoscalar (VPP) vertices appearing in the diagrams 
in Fig. 4 and then predict their couplings in the FMV model. 

The most general 0{p^) weak octet VP'y vertex for the processes of our interest has 



been already worked out in Ref. |T^ and it is 



with 



and 



CwiVPj) = GsF^iV^j;^), (39) 



= E , (40) 



rr^ = {«^A/f}, 

= {{A,n^},/f}, 
Tg""^ = {Au^ff). (41) 



In the same way the most general 0{p^) weak octet VPP vertex, without including mass 
terms, is 

CwiVPP) = GsF^iV^'^j;^,), (42) 

with 

J!^u = E , (43) 

and 

Si, = iu^Au,, (44) 
= ^t^Uv ) ■ 

In Eqs. (^OD and (^), the Kj and cTj are the a priori unknown coupling constants to be 
predicted in a model or given by the phenomenology when available. 

With the weak vertices in CwiVP--^) and CwiVPP) and the strong/electromagnetic 
vertices given by Cv in Eq. (0) we can compute now the diagrams in Fig. 4 and we get 
the direct contribution to the magnetic amplitude as given by vector meson dominance 
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7 

6) 




Figure 4: Diagrams contributing to the direct vector meson exchange magnetic amphtude 
M^l The empty box is generated by Cw{VP-r) in Eq. (H) or Cw{VPP) in Eq. (H) 
and the black circles by the strong/electromagnetic lagrangians £y in Eq. (|13). In b) and 
d) the crossed vr^ ^ vr^ diagrams are also understood. 
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in terms of the and cij couplings 

(6) _ 16 Gsem% 



dir 



16 Lr^errij^ ( 
—= — 1 9v [S/ta + 3^3 - (5/t2 + <0K:i) z^. 
3v2 'myt'T, L 

-\/2hv [2cri - (ScTi + 0-2) -^3 



(45) 



In Ref. [0 we have proposed a Factorization Model in the Vector couplings that has 
proven to be very efficient in the understanding of the vector meson contributions to the 
K — > 7r77 and Kl —>■ 77* processes. The prediction for Ki in the FMV model has been 
worked out in Ref. [0]. To evaluate the CTj we proceed analogously by applying factoriza- 
tion to construct the weak VPP vertex. A detailed evaluation is shown in Appendix B. 
The predictions for the coupling constants in this model are 



^FMV 







FMV 



K 



FMV 



{2hv - iv) VVP-y 



VvPf 



a 



FMV 
1 



fv [ , AT 

-j= \ r]v + 4:Lgj^7]vpp 



„FMV 
"^2 



FMV 



FMV 



2iv Vvp^ 
2iv r]vp-y 



FMV 




(46) 



my ^ 
VV + 2L9-^?7ypp 



where 



16y27r' 



ml 
:JvY2 



(47) 



and the last identity is exact in the Hidden gauge model and also well supported phe- 
nomenologically. As explained in Ref. [Q, iy is a contribution given by the O(p^) WZW 
anomaly Eq.(|^. In Eq. (^) the unknown factorization factors are rivp-y that comes from 
the 0{p^) weak VP'j vertex, rjvpp from the 0{p^) weak VPP vertex and rjv from an 0{p) 
weak PV vertex (see also Appendix B) and that we already have met in the predic- 
tions for NY^^ in Eq. (|25|). Naive factorization, however, would put rjvp-y = rjvpp = rjv- 
In Ref. [14 1 we have shown that the phenomenology of — 7177 and Kl 77* is well 
described with rjyp^ ~ 0.21 as predicted by the naive Wilson coefficient in Eq. (^) (i.e. 
no enhancement of the AI = 1/2 transitions) and we will use this result in our numerical 
study. For rjypp and rjv we do not have still any information. However we will assume 
rjvpp — Vvp-y because of factorization and we will leave rjv free. This is because while 
r]vpp and r]vp-y appear both at 0{p^) in the weak chiral lagrangian, rjv appears at 0{p) 
and therefore could differ appreciably from the naive result. 

With the FMV predictions for the couplings Ki and ai we can give the direct M^l 
amplitude in this model. Substituting in Eq. (^) we get 

(6) _ Ggem^ 



M, 



dir,FMV 



■rv 



(Wp7 + Vv) (1 - 0^3) 
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+ WP7 7T^ (2-5^3) (4^ 
4:hv 



+ 2rivpp Lg-^ (2 - 5^3) 



where ry has been defined in Eq. (|37D, iy in Eq. p7| ) and we have put = 0. Note that 
M^l gets a dependence on the photon energy. 

Another procedure to implement the hypothesis of factorization has also been used 
in the literature 0, In this case one determines first the strong action generated by 
vector meson dominance and then applies the FM as given by Eq. (^). This is the so 
called procedure A [0, ^ in opposition to the FMV model where one first determines in 
Factorization the weak couplings of vectors and then integrate them out (procedure B). 

According to the method A the vector generated 0{p^) strong action necessary to 
compute the direct magnetic amplitude for Kl 7r+7r^7 has two contributions : i) 
the one given by integrating the vector mesons between the terms in gy and hy in Cy 
(Eq. (|10D), and ii) a term obtained by integrating the vectors between the terms in fy 
and hy in Eq. (|l^). Once factorization is applied these generate the following magnetic 
amplitude : 



(49) 



(2-3^3) + ^(2-5^3) 

^gv 



where ry has been defined in Eq. (|37|) and kp is the factorization parameter. In Mf^^pj^ 
the term proportional to fy is the one given by the second contribution (ii) above. In 
Ref. only the term generated by the first structure (i) has been considered. 

The theoretical comparison between the two models motivates the following consider- 
ations. 



1) As has been shown in Ref. the FM and FMV models do not give the same 
contributions. In fact the FMV model recovers the results of the FM but it is able 
to generate a more complete set of effective operators. The fact that one recovers 



more chiral operators in the FMV model was shown to be important in Ref. [14| for 
the understanding of the vector meson contributions to K 7177 and Kl 77* . 
In Kl —* 7r+7r^7 the weak VPP vertex generates an extra contribution (last line in 
Eq. (HD) that the FM does not get [0] 0. 



2) The af^^ couplings in the weak VPP vertex have two different kinds of contribu- 
tions. From Eq. (^6|) we see that there is a term proportional to rjypp and a piece 



^In comparing our result M^^^ fmv with the one given by the FM in Ref. § care has to be taken 
because the anomalous contribution given by the term proportional to Lg in the FM (Eq. (5.24) of 
Ref. 0) corresponds to our term proportional to £y in Eq. ( |4^ ) and coincides with it once the relation 
L9 = fvgv/'^ = /v-/4 is used. 
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proportional to rjy. The first one, analogous to the 'qvp^y contribution in nf^^, 
comes from the application of factorization through the left-handed currents in the 
direct VPP vertex. As explained in Appendix B the origin of the term proportional 
to riv is a weak shift in the kinetic term of the vector fields in the conventional 
vector formulation. This framework, that we have proposed in Ref. [^, has been 
able to generate new features in the 0{p^) Ni couplings of £4^*^'"^ like, for example, 
non-vanishing spin-1 resonance exchange contributions to the N28, N31 couplings 
in Eqs. ( p!8| , p5D that are predicted to be zero in the FM with the antisymmetric for- 



mulation [|I^ . It is reassuring to notice that the term proportional to r/y allows us to 
recover the result of the FM |0] for 7r+7r"7 (that is 0{p^)) in agreement with 

our claim that the procedure we are applying uncovers new structures but does not 



miss any contained in the previous FM. This shows that our approach in Ref. 
is also consistent at 0{p^) and hence with the chiral expansion. 

We think that the FMV model seems to give a more complete description and therefore 
we will analyse the Kl — > 7r+7r~7 process in this model. 

Using the relations iy = 4:hv (Eq. (|47|)) and Lg = /y/4 we can give a simplified 
expression for M^^^j^y, 



,.(6) Gsemj^ 

2-K 



2u + u' ~ ^(13 + 10cu + 3cj') Z3 



(50) 



where uj = rjvpp/'ijvP'y, ^' = Vv/vvP'y and naive factorization would suggest u = u' = 1. 
However, as commented before, we would like to emphasize that while the value of rjyp^ ~ 



0.21 has been fixed in Ref. [14] from the phenomenology of K ^ Tcyy and Ki — > 77* , no 
similar constraint affects rjypp or t]v and it could happen that, in contradistinction with 
rjvPf, these weak vertices are affected by the AI = 1/2 enhancement and therefore u and 
u' could differ appreciably from the unity. It is worth to notice that, in principle, both 
contributions to M^^^_^mv coming from the weak VP'y (the pure numerical term in the 
brackets in Eq. (|50D) and VPP vertices (terms in uj and u') are quantitatively comparable 
and of the same order. Thus we find that the phenomenological assumption by Heiliger 
and Sehgal |^6|] that the vector exchange generated DE amplitude of Kl —>■ 7r+7r^7 is 
dominated by the weak VP'y vertex as in i^' — * 7177 and — > 77* , is not supported by 
our analysis. 

Moreover, comparing M^^^ pj^y and M^J^ in Eq. (^) we also see that both direct and 
indirect contributions are comparable as it happens in K Tc-f-f and Kp 77* . 

4.2 Loop amplitudes 

The leading loop contributions to the magnetic amplitude of Kp — > 7r+7r^7 are 0{p^) 
in xPT . There are two different kinds of loop generated amplitudes. 
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The first type corresponds to the Feynman diagrams in Fig. 5 which contribution is 
given by the lagrangian densities c''^^^~^ in Eq. (|T^), the strong/electromagnetic £2 in 
Eq. and the WZW anomalous action Zan[U, i, r]wzw in Eq. (|^). Thus the main feature 
is that this contribution is completely specified on chiral symmetry and anomaly grounds 
and, therefore, it is model independent. 

In Fig. 5. a the fields running in the loop are the pairs {rj^^n^) and (tt*^,??^). In the 
pole-like diagrams in Figs. 5.b-5.e we keep those where the external kaon leg undergoes 
the weak transition to 7r°, r^g and disregard those (suppressed by m\/m\) where the 
transition happens in a final leg. In the diagram in Fig. 5.b only is running in the 
loop and its final contribution vanishes because a cancellation between the 7r° and the r^g; 
in Fig. 5.C we have tt^,K^,K^\ in Fig. 5.d the couple (7r^,7r°) and in Fig. 5.e we have 
vr^, vr'^, and . We note also that there is no pole contribution with r^g in the diagram 
in Fig. 5.d. In Fig. 5.f-5.h the weak transition happens inside the loop. In the diagram 
in Fig. 5.f we have (1,2,3) = (Kg, 7r°, i^'i), {Ks,ris, Kl) ; in Fig. 5.g we have (clockwise 
running in the loop) (1, 2, 3) = {K^, n~,r]s), {K^ , vr", 7r°), {Kl, 7r°, 7r+), {Kl, ?7g, 7r+) and, 
finally, in Fig. 5.h the {K'^,tt'^) and {K^,7i^) pairs appear. 

The result is divergent and in the MS subtraction scheme is [| 

,^(6) _ C^serrij^ f rriK 



9 - 9 3(1 -r^) V/^V 3(1 -r^) \ ^^ 



+ 



"5 














-> 




3 











Z3 



2(1 -r^ 



+ K[{p - p+f ,m\,ml] + K[{p-p^f,m\,ml 

+ K[{p+ + p.)\ m% ml] + K[{p+ + pJ)\ ml, 

1 - 1 - 

4 4 



2 

+ -rl [F[{p-p+)\ml] + F[{p-p 



where the functions F[g^,m^] and K[q'^,ml,m2] have been defined in the Appendix C. 

We have compared the loop evaluation of 7r° 7r+7r~7 and ?7g 7r+7r~7 (necessary 
for the pole-like diagrams in Fig. 5) with the results shown in the Ref. ||3^ and found 
complete agreement. From our Eq. ( ^l]) it is easier to obtain the appropriate momentum 



dependence that is not properly specified in Ref. [37|. Moreover the divergent part of 



^We recall that we have neglected diagrams with a weak transition in a external pion leg (suppressed 
by r^) and these cancel the (1 — r^) in the denominators of M^^^yy. 
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9) h) 



Figure 5: Diagrams contributing to the loop magnetic amplitude Mp^^^y. The empty box 
is generated by C^^^^~^ in Eq. (|l^ , the white circles by the WZW anomalous action Zan 
in Eq. (|^) and the black circles by the strong/electromagnetic lagrangian C2 in Eq. (^. 
In a), d) and g) the crossed 7r+ ^ tt^ diagrams are also taken into account. The particle 
content in the loops is explained in the text. 
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these diagrams comes from the strong/electromagnetic loop because the weak transition 
happens in the external leg. We have checked, then, that the strong intrinsic parity 
violating 0{p^) lagrangian can absorb these divergences. 

The second type of loop contributions corresponds to the Feynman diagrams in Fig. 6. 
Their basic characteristic is that the weak vertices are generated by £4^'^'"^ in Eq. (IT 



and therefore the generated amplitudes are proportional to the, a priori unknown, Ni 
couplings. The strong vertices in Fig. 6 are generated by £2 in Eq- (ID- In Fig. 6. a 
the fields running in the loop are the pairs (7r+,7r~) and {K~^,K~)] in Fig. 6.b we have 
7r+, TT^, ris, and finally in the graph in Fig. 6.c there are the pairs tt^), (-ft'^, vr°) 
and {K"^,!]^). We also have taken into account the field and renormalization due to 
the leading M'^^^ amplitude in Eq. 



As in the case above the result is divergent and using the MS scheme it reads 



(6) 



3 In 



m 



K 



+ (3iV29 - N- 



30 J 



+ 2iN29 + Nsi) 



2 In 

+ K[ip 
+ K[ip- 
5 

6 + 



m 



K 



/i 



P-] 



+ K[{p - p^ 



K[{p-p+f,m 



P- 



-r^ + - ) In ( -f 1 + - In 




In 




4 2 

-rlF[{p++p_Y,ml] - 



m 



K[{p - p+f ,m\,ml] +K[{p-p. 



\2 2 2l 



(52) 



^3 



and the functions F[q^^m?] and K[(f 
r(6) 



, , rn^ have been defined in the Appendix C. 



We notice that M\j', depends on the four A^28) ^31 unknown couplings in the 0{p 
weak lagrangian, but only three combinations of them, A^28, 3A^29 — A'so and A^29 + A'ai, 
appear. 

In our numerical study of the branching ratio and spectrum in the next section we 



have worked with the full expressions of M^^\w and M);/ with /i = nip. However in 
order to make more transparent our study of the slope in Kl tt^ti^^ (Eq. (0)) we 
have linearized the real part of the loop contributions by Taylor expanding the amplitude 



r{6) 
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T + 




Figure 6: Diagrams contributing to the loop magnetic amplitude Mj^r, . The crossed box 
is generated by £4^*^'"^ in Eq. ([I8[) , and the black circles by the strong lagrangian C2 in 
Eq. d^). In c) the crossed tt^ ^ tt~ diagram is also considered. The particle content in 
the loops is explained in the text. 
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assuming 4r^ < -23 <^ 1- We get for the full 0{p^) loop contribution 

(6) ^ Gsem% / rriK 

Loop - ^^2p^ \^^p^ 



MUop = ^ 2f UriT [fli + ^L^s] , (53) 



where 

= 0.54 - 2.63 ai - 2.07 (02 + 204) - 0.13 (02 - 203) 
Bl = 0.71 + 0.67 (02 + 204) - 0.07 (02 - 203) . 



(54) 



The linearizations give values for Re{M^^^p) in very good agreement with the exact results 
(in a very few percent inside the domain of 2:3). From Eq. ( |5^ we notice that if ~ 0{1) 
(as expected) the main contributions come from the terms in ai oc N2S and (02 + 204) oc 
{N29 + N31). 



5 Analysis of Kl n^n 7 



Collecting all our previous results for the magnetic amplitude of 7r"'"7r~7 (with 

the linearized version of the loop contributions) we get for the factor m and the slope c 
defined in Eq. (^2]) the results 



m 



02 + 2a4 - Fi + ry [ 1 + rjvpj (3 + 2 + 



UJ 



dL 



1 

m 



3 + 



VVPy 



■l3 + 10cu + 3cj') + 



(55) 
(56) 



We recall here what is known and what is not in the expressions above: 



a) The O(p^) 02 + 204 term is unknown. In Subsection 2.2 we have given the predictions 
for Oj in the FM (anomaly) and the vector formulation framework using factorization 
(spin-1 resonances) in terms of one free parameter 77. We will study our observables 
in a reasonable range of this term : < 02 + 204 < 3 (we remind that factorization 
predicts Oj ~ 0{1) and positive). 

b) -Fi has been defined in Eq. (^) and the only parameter we will allow to be free is 
the one measuring the breaking of the nonet symmetry in the weak vertex p defined 
in Eq. (jl^). We will consider < p < 1 as a conservative working region. 

c) ry is given in Eq. (P?]); rjvp-y has been fixed in Ref. |T^ to rjvp-^ — 0.21 as commented 
before. 
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/7o -1- '\)n A 


n 2 


4 


fi 


8 


1 


1.2 




— U. iO 


— U.oO 


— i.OO 


9 9/1 


9 Q/1 


— o.oo 


bL 


0.85 


0.98 


1.12 


1.26 


1.40 


1.54 



Table 1: Loop contributions ai and 6/, in Eqs. ( p3| , |5^ for different values of the a2 + 2a4 
parameter as explained in the text. 



d) The ratio u = rjypp/rjvp'y, defined in Eq. (|50|) , is also free and we will allow a 
domain 0.5 < u < 1.5 coherently with the naive estimate of factorization that sets 

~ 1. However we will keep free uj' = r)y/rjyp^ and will allow a reasonable range 
of values for r/y, i.e. < r^y < 1. 

e) The and 6^ parameters have been defined in Eq. (Q) and therefore they depend 
on the unknown couplings. Our procedure will be the following. For every 



value of the 0{p^) a2 + 2a4 free parameter we use the results in Eqs. (^ , |25|) with 
Vv = Van = ?7 to fix a value for rj as given by 

1 = . (57) 

3 + —^fyhy 

With this value then we fix a\ and 02 — 2a3 that appear in the loop contributions. 
For definiteness we have not included the axial contribution in A^j and therefore in 
Eq. (|57D. If included it would add constructively in the denominator giving a smaller 
value of 77. From Eq. (|54D we see that the dependence of the loop contribution in 
02 — 2a3 is very mild and only a\ is relevant. Values for ai^ and 6^ for fixed + 2a4 
are collected in Table 1. We notice that the dependence on the combination 02 + 204 
is bigger in a^. 

We have therefore three, a priori, free parameters : p, lo and 02 + 204. Using now 
the experimental values for the branching ratio and the slope of Ki^ -n^-n^^ we can 
proceed to the analysis. However there is another process where the breaking of the nonet 
symmetry in the weak sector could play a role : — 77 . A short note on this process 
is now needed. 
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Ki — >■ 77 is a decay dominated by long-distance contributions which first non vanish- 
ing amphtude is 0{p^) in xPT P^ , |D|. Neglecting CP violating effects the amplitude is 
defined by 



A [KLip) 7(^1, £1)7(92, 62)] 



i A^^ e^^^r ei (gi) £2(^2) Qi qI 



(58) 



There is, up to now, no full computation at 0{p^) in xPT . At present the amplitude is 
taken 



A. 



77 
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Gs F-n OLem -^2 



(59) 



where 
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3(1 



[(1 +Ocos^ + 2V2psin^] 



3(1 -r^o 



[2V2pcos0- (1+Osin^] 



^ sin Q 



^co^Q -1^P&- sin^ 
F 

2^2-^ cos e 

Fq 



(60) 



comes from the pole-like diagrams with the same definitions than Fi in Eq. ( pSf ) [3. In 
our previous study on Kl —>■ 77* we noticed that, keeping only the term in F2 in 
Eq. ( |55| ) with p = 1 and F^^ = F^ = Fq all the vector meson dominance models known 
were giving a wrong sign for the slope in Kl 77* . Due to the poor knowledge of 
the Kl — > 77 amplitude we think that the problem could be the sign of this amplitude. 
Using the experimental figure \A^^^\ = (3.51 ± 0.05) x 10~^GeV~^, from Eq. 



we get 1-^2^^! = 0.883 ± 0.013. We have fixed the range < p < 1 by the conservative 
assumption that iFg'^^'"'! does not differ from the experimental figure more than a factor 
2. We will include this process in our study and will combine, therefore, 77 and 

Kl — > 7r+7r^7 . 

In our numerical analysis we use the following values for the known parameters: the 
mixing angle between 773 and r]o that appears in Fi and F2 is taken 6 ~ —20°; the SU{3) 
breaking parameter ^ defined in Eq. (^) is ^ ~ 0.17 [|^; the ratios between the decay 
constants are F^/Fg ~ 0.77 and F^/Fq ~ 0.98 || □ the couplings defined in Eq. 
are fv = 2gv ^ 0.18, hy ~ 0.037 and By = 2/iy. 

Our analysis shows the following features : 

1/ There is a very mild dependence in the u parameter for the indicated domain 0.5 < 
uj < 1.5. It gives an uncertainty in our predictions on Kl 7r+7r~7 of around 
10% for the branching ratio and 5% for the slope. Hereafter we take = 1 in our 
analysis. 



loin Ref. 



evaluation using £4'^'^' ^ 



1[ has been pointed out an extra contribution to A.yj coming precisely from a one loop 
We have not been able to recover that result and found the claimed new 
contribution vanishes, therefore we do not take it into account. After reviewing their calculation the 
authors of Ref. |^ now agree with us. 

These values are consistent with the recent analysis of Ref. p3]. 
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2/ In Table 2 we show, for fixed values of 02+204, the region in p that gives a reasonable 
result for B{Kl 7r+7r^7 ]E* > 20MeV)£)E, the prediction for the slope c in 
Ki 7r^7r^7 and the values of F2 in A^^. We see that the branching ratio can 
always be accommodated for reasonable values of 02+204 and p. Taking into account 
IF2I in Kl — > 77 , however, there is a restriction on the lower and upper limits of 
02 + 204, roughly 0.2 < 02 + 204 < 1.2 and also an internal region of this interval is 
excluded because IF2I is too small. Thus the combined analysis of the experimental 
widths of Kl —>■ 7r+7r^7 and Kl — > 77 are consistent in two regions of parameters 
that are ruled only by the value of the combination 02 + 204 : 

Region I : 0.2 < 02 + 204 < 0.6 , < p < 0.4 , 
0.05 < ri < 0.13 , F2 < ; 

Region II : 0.8 < 02 + 204 < 1.2 , 0.6 < p < 1 , 
0.18 < T] < 0.23 , F2 > . 

We remind that pion pole dominance in Kl —>■ 77 demands F2 > 0. Our analysis, 
however, shows that both possibilities are consistent in the framework we have 
developed here, and that the breaking of the nonet symmetry in the O(p^) weak 
vertex could consistently be as big as p ~ 0.2. Notice that if we had included the 
axial-vector contribution in the determination of 77 in Eq. (^) this parameter would 
diminish by a 15%. 

3/ As can be seen in Table 2 the predicted slope c of Kl —* 7r"'"7r~7 is in a small 
range c~ —1.6, —1.8, independently from the range of parameters and in excellent 
agreement with the experimental figure Cexp = —1.7 ± 0.5. As seen from Eq. ( p^ 
the slope has two contributions : vector exchange and loops. The role of the latter 
is crucial in stabilizing the prediction. This can be seen in Table 3 where we show 
the values of the slope including only the vector amplitude (c ~ —1.5, —2.8) and 
the total result obtained by adding the loop contribution (c ~ —1.6, —1.8). The 
corrections to the width coming from the loop amplitudes are mainly relevant in 
Region II where can reach even 100% for 02 + 204 = 1.2. This is because in this 
region there is an almost complete cancellation of the term 02 + 204 — Fi of m in 
Eq. (|55|) and thus providing instability to the result; axial- vector exchange or higher 
0{p^) corrections, not included in this analysis, might then be relevant. In Region 
I the loop contribution to the branching ratio is never bigger than 20%. 

In Fig. 7 we plot the spectrum of the direct emission Kl 7[~^tt~^ in the photon 
energy (E*) for the two extreme values 02 + 204 = 0.2 (Region I) with p = 0, c = —1.61, 
and 02 + 204 = 1.2 (Region II) with p = 1, c = —1.78. The tiny difference between both 



26 



a2 + 2a4 


P 


B{Kl 7r+7r-7 ; E* > 20MeY)DE x 10^ 


c 


F2 


0.2 


-0.05 


4.09 


-1.52 


-1.89 


0.00 


3.19 


-1.61 


-1.72 


0.05 


2.93 


-1.71 


-1.55 


0.4 


0.15 


4.00 


-1.56 


-1.22 


0.20 


3.40 


-1.64 


-1.05 


0.25 


2.85 


-1.71 


-0.89 


0.6 


0.35 


3.93 


-1.59 


-0.55 


0.40 


3.34 


-1.68 


-0.39 


0.15 


2.80 


-1.78 


-0.22 


0.8 


0.55 


3.85 


-1.62 


0.11 


0.60 


3.26 


-1.71 


0.28 


0.65 


2.72 


-1.82 


0.45 


1.0 


0.75 


3.78 


-1.65 


0.79 


0.80 


3.20 


-1.75 


0.95 


0.85 


2.67 


-1.85 


1.11 


1.2 


0.95 


3.70 


-1.68 


1.45 


1.00 


3.12 


-1.78 


1.61 


1.05 


2.60 


-1.89 


1.78 



Tabic 2: Range of values of p that give a reasonable prediction for the branching ratio 
and the slope of 7r+7r^7 at fixed value of 02 + 204 for = 1. The prediction 

for F2 is also given. The experimental values are B{Kl n^n^'j ; E* > 20MeV)DE = 
(3.19 ± 0.16) X 10-^ Cexp = -1.7 ± 0.5 and IF2I = 0.883 ± 0.013. 
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n 2 


4 


n 6 


n 8 


1 n 


1.2 


P 


n n 
U.U 


n o 
U.z 


n A 
U.4 


n R 

U.D 


U.o 


1 n 
i.U 


c [vntioi b oiiiy ) 




— i. ( u 


1 SQ 
— i.oy 


— z. iO 


9 /I 


9 89 
— z.oz 


c (vectors + loops) 


-1.61 


-1.64 


-1.68 


-1.71 


-1.75 


-1.78 



Table 3: Comparison between the values of the slope c of the magnetic amplitude of 
Kl 7r+7r^7 for a fixed value of 02 + 204 as given by vector exchange only and adding 
the loop contributions in the total result. When the latter are included the prediction is 
much more steady. 



curves shows the model dependence of our prediction. We have normalized the spectrum 



to the 1937 events reported by the E731 experiment in Ref. ||3§| that are also displayed. 
From this Fig. 7 it would seem that a bigger value of the magnitude of the slope (and thus 
shifting the spectrum to lower values of E*) might be needed. We would like to stress 
that the apparent disagreement between the data and the theoretical predictions is due 
to the fact that data seem to support a larger slope. However the form factor fitted by 
the experimentalists gives a slope in very good agreement with our prediction. A more 
accurate experimental determination is required (the events at low energy seem to have 
a huge error and thus a better background subtraction from bremsstrahlung is needed). 

In Ref. a parameter / measuring the fraction of the AI = 1/2 K ^ vrvr amplitude 
due to penguin effects was introduced. This ratio can be related with the parameter p 
introduced in Eq. (^) that measures the breaking of the nonet symmetry in the weak 
sector through the relation (assuming CP conserved) p = (3/ — l)/2 |^]. There are a lot 
of uncertainties in the determination of / and little can be said. We just recall that for 
the generous range 0.1 < / < 0.8 it implies —0.2 < p < 0.7. Coming back to our analysis 
we see that Region I implies 0.3 < / < 0.6 while Region II allows 0.7 < / < 1 in the 
boundary of the reasonable range. 

A joint analysis of Kl — > '7r+7r~7 and Kl 77 has also been carried out previously 
by other authors. Cheng ^ only includes the term 02 + 2a4 — Fi in m (Eq. (|55|)) and, 
moreover, arbitrarily input 02 = 04 = 1. His conclusions are that the DE — > 7r+7r~7 
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Figure 7: Spectrum in E* {E* > 20 MeV) of the Kl 7r+7r^7 (DE) process normalized 
to the 1937 events of the E731 experiment of Ref. [BB] that are also displayed. The 
continuous line corresponds to 02 + 2a4 = 0.2, p = and u = 1 with a predicted slope 
of c = —1.61. The dotted line corresponds to 02 + 2a4 = 1.2, p = 1 and a; = 1 with a 
predicted slope of c = —1.78. 
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branching ratio excludes a small value for the p parameter. Ko and Truong |46] and 
Picciotto have considered the contributions of the reducible anomalous amplitude 
(term in Fi) and the indirect Yector interchange (keeping, however, the whole vector pole 
structure), but not the leading 0{p'^) amplitude : they assume 02 = 04 = 0. These 
authors conclude that only a small breaking of the nonet symmetry in the weak sector is 
compatible with their analysis (p ~ 0.8, 0.9) and therefore, that the vr'' pole dominates 
the Kl — > 77 amplitude. Moreover their results give a small dependence of the magnetic 
amplitude in the photon energy and, consequently, a too small slope c. 

All these previous analyses emphasize that only our Region II is allowed. In our more 
complete analysis we show that this is not the case and both solutions are kept open. 
In fact several features (the sign of the slope in — > 77* and the ratio /) lead to the 
scenario described by the parameters in Region I. 

A comment about the convergence of the chiral expansion is also interesting. The bulk 
of the contribution to the width of Kl 7r+7r^7 comes from local contributions at 0{p'^) 
and 0{p^) , but the loop amplitudes are relevant for the region 0.8 < 02 + 204 < 1.2. The 
slope in E* starts at 0{p^) in xPT and it is mostly due to the vector exchange but the loop 
contribution can be as big as 30% and in any case is crucial to give a steady prediction. 
Local 0{p^) corrections due to resonance exchange (other than vectors), even if suppressed 
by heavier masses, could give a non-negligible contribution to the width if numerical 
factors overcome the suppression. The analysis of 0{p^) is very much cumbersome. It 
is curious to notice, however, that the 0{p^) (and higher orders) amplitude specified in 
Eq. (|38|) (that we have not included in our analysis) does not seem smaller than the 
pure 0{p^) M-^\ in Eq. (|37| ) because, in principle, Fi ~ 0{1). If the chiral expansion 



has to converge one would expect either a cancellation with other 0{p^) contributions or 
that the parameters allow a self-suppression. Hence for our preferred values of 02 + 204, 
p = 0.8 implies Fi ~ 0.96 but p = 0.2 gives Fi ^ —0.07 that is one order of magnitude 
smaller. This fact together with the statement pointed out above that in Region II the 
0{p^) loop corrections to the width are huge indicates that a better convergence of the 
chiral expansion (in the terms considered here) seems to be achieved with 02 + 2a4 and p 
small [5 



^^We have studied the scale dependence of the loop amplitude. As a matter of fact this is not negligible 
by itself (i.e. ^ 30% in the loop for a reasonable range of variation of /x). However the loop contributions 
play a very marginal role in the observables (i.e. rate and spectrum) for Region I; thus this change is 
negligible. This is not the case in Region II (where strong cancellations make the final result sensitive to 
any small change) confirming our conclusion that predictions in this region are rather loose. 
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6 Conclusion 



The measured energy dependence of the spectrum for the process Kl —>■ 7r^7r^7 shows 
that the slope of the magnetic amphtude defined in Eq. ( p2D has a relative big size thought 
to be due to a relevant vector meson exchange contribution at 0{p^) in xPT . Moreover 
with our present knowledge of the phenomenology the width of the process cannot be 
predicted in a model independent way. In Ref. we pointed out that the sign of the 



Kl 77 amplitude has to be opposite to the one predicted by pion pole dominance 
in this decay. Both processes are correlated because receive a reducible contribution 
depending on the unknown p parameter (Eq. (|1^) that measures the breaking of the 
nonet symmetry in the 0{p'^) weak vertex. 

In this work we have presented a systematic study of the magnetic amplitude of the 
process Kl 7r+7r^7 up to in xPT . We have computed the local contributions 

(due to a reducible anomalous amplitude and vector meson exchange) and the one-loop 
magnetic amplitude at this order. 

The 0{p'^) magnetic amplitude (M^'*-') and part of the loop contribution (M^.^) are 
model dependent because the couplings Ni in £4^'^'"^ Eq. ( p!8D are not known from the 
phenomenology. We have included the spin-1 resonance exchange contribution to the 



relevant Ni couplings that we have evaluated in the approach described in Ref. |]26 
Moreover part of the vector exchange amplitudes due to diagrams in Fig. 4 where a weak 
I/P7 or VPP appear are also not known. Driven by the good description we have achieved 



in Ref. |14] for the processes K —> 7r77 and Kl — > •yi^i we have used the Factorization 



Model in the Vector Couplings (FMV) [Tj, ^ in order to provide the direct weak VP'y 
and VPP vertices. 

The conclusions of our analysis can be stated as follows : 

1 / We are able to accommodate the experimentally determined branching ratio Bt^Kl — 
7r^7r~7 ; E* > 20MeV)DE inside a reasonable range of variation of the parameters 
02 + 2a4 and p. When combined with the Kl — > 77 amplitude we find two regions 
for the combination 02 + 204 : 0.2 < 02 + 204 < 0.6 and 0.8 < O2 + 204 < 1.2. Away of 
this range l^-y^l is predicted to differ for more than a factor two of the experimental 
figure. Remarkably in any of these cases the value given by Eq. (|5^) is 77 < 0.3 im- 
plying, therefore, that the octet operators are not enhanced and are consistent with 
the perturbative expectation provided by the Wilson coefficient [rj ~ 0.2). This we 
can interpret as a fiducial behaviour of factorization showing that the dynamical 
features introduced in our analysis through this hypothesis are compatible with the 
perturbative expansion that becomes predictive a this point. 

We have also concluded that the Region II is perhaps sensitive to effects not in- 
cluded here like 0{p^) axial- vector exchange or non-resonant contributions and 
0{p^) corrections. 
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2/ Independently of the parameter variation we consistently predict a small range for 
the slope of Ki —>■ 7r"^7r~7 , roughly c ^ —1.6, —1.8 in very good agreement with 
the experimental result Cexp = —1.7 ± 0.5. The loop contributions are shown to be 
crucial to stabilize the slope in this small range. As a consequence the spectrum 
plotted in Fig. 7 is a pure prediction of our framework. Corrections to it can come 
from other resonance interchange at 0{p^) (axial-vectors, scalars, ...) or 0{p^) 
effects. 

3/ From Table 2 we see that, in our joint analysis, both signs of the 77 ampli- 

tude are consistent with the phenomenology, in contradistinction to the conclusion 
reached by previous studies. Hence the result F2 < (pion pole dominance implies 
F2 > 0) is allowed in consistency with our statement in Ref. |T^ that the slope in 
Kl 77* , experimentally determined, would imply that sign. If this is the case we 



would conclude that, in opposition to previous results [§3, ^ a bigger breaking 
of the nonet symmetry in the weak vertex is called for (p ~ 0.2). Of course, at this 
point, it cannot be excluded a solution with p ~ 0.8, and blame the change of sign 
of the Kl — i> 77 amplitude to a, still unknown, huge higher order correction. 

The relevance of other resonance exchange generated 0{p^) or higher order contribu- 
tions would be very much clarified once a more accurate measurement of the spectrum in 
E* is carried out. The experiment KLOE at DA$NE expects to have 35, 000 — 7r+7r^7 
events p/year due to DE with E* > 20MeV and therefore should be able to improve 
statistics and accuracy in this channel. Other two new experiments NA48 at CERN and 
KTEV at Fermilab could also bring new features on this process. These are going to focus 
on the related 7r+7r^7* process where to keep under theoretical control the leading 

real photon channel is crucial to disentangle its dependence. 
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Appendix A : Kinematics on Ki tt^tt 7 

In Section 3 we defined the amplitudes and wrote explicitly the differential cross 
section for an unpolarized photon in Kl n^n^'j . Here, for completeness, we collect 



some useful kinematical relations Iq, 33 



In Kl 7r+7r^7 the most useful variables are : i) the photon energy in the kaon rest 
frame {E*), ii) the angle (6) between the 7 and vr"*" momenta in the di-pion rest frame. 
The relations between {E*, 6) and the Zi defined in Eq. (|28|) are 

^± = 77^(1 T/3cos^), z, = ^, (A.l) 

IrriK rriK 



where (3 = Jl — 4:711^ /{m\ — 2mKEZ). The kinematical limits on E* and 9 are given by 



< < ^""^ , -1 < COS0 < 1 . (A.2) 



The differential rate in terms of these variables for an unpolarized photon is 

a'r {Elf (3^ ( 2e:; 



dE*dcose 5127r^mi. V 



1 1 sin^^ f \E\^ + IMP) . (A.3) 



Appendix B : The weak VPP vertex in the Factoriza- 
tion Model in the Vector couplings 



Following the same procedure that we used in Ref. |jTj] to construct the weak VP^y 
vertex we proceed to evaluate the weak VPP vertex for — > 7r+7r~7 . 

The bosonization of the Q- operator in Eq. (^) can be carried out in the FMV from 
the strong effective action S* of a chiral gauge theory. For later use let us split the strong 
action and the left-handed current into two pieces : S = Si + S2 and J'^ = + JT"^, 
respectively. Then in the factorization approach the Q_ operator is represented by 

Q- - 4 [{x{ji,jn) - {^Ji){j2n - {^j'){jr)] , (b.i) 

with A = (Ae — iXr)/^ and, for generality, the currents have been supposed to have non- 
zero trace. 

In order to apply this procedure to construct the factorizable contribution to the weak 
0{p^) VPP vertex we have to identify in the full strong action the pieces (at this chiral 
order) that can contribute. It turns out that there are four terms in the strong action 
that can play a role. Analogously to the specified procedure we define correspondingly : 



5 



S — Sy + Sp 

Sv = Sy-y + Sypp , (B.2) 

c _ CX I CX 

Dp — D2 + O4 , 
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where the notation is self-exphcative and the actions correspond to the lagrangian densities 
proportional to fy and gy in Cy Eq. (jIDD, £2 in Eq. and £4 in Eq. (|^). 

Evaluating the left-handed currents and keeping only the terms of interest we get 

6Sv fv 



6F 
6Sp 



V2 



:b.3) 



2 

Then the effective action in the factorizable approach is 



4 Gg Vvpp 



j6Sv ^Sp\ 6Sv,,SSp 

^M^'^r ~ ^ ^^^^^ 



+ h.c. 



(B.4) 



There is however an extra contribution that remains to be taken into account. In Ref. 
we have shown that in the conventional vector formulation one can construct a weak 0{p) 
coupling involving vectors as 

'ujX{^{V„u^}) +uU^u,){Vn] , (B.5) 



where the couplings ujy are in the FM 

2 



UJ 



V 



V 2 — ^ fv Vv 



V 



V 

1 > 



:b.6) 



with r]v the unknown factorization factor. If one rotates away these terms through the 
shift 



V, - 



GgFt 



LuX {A,M^} + uX (An^) 



the kinetic term of the vector mesons 

1 



C 



m 



K 



(B.7) 



(B.8) 



V^y has been defined in con- 



generates also an 0{p^) weak VPP coupling. In Eq. 
nection with Eq. ([lO|). 

Comparing the full lagrangian density Cv/{^PP) in Eqs. we find the pre- 

dictions of the FMV model for the couplings cTj. Using fv = 2gv we get 



FMV 
1 



V2 



r]v + 4 Lg r]vpp 



FMV 
^2 



-V2f 



vVv 



:b.9) 



FMV 



V2fv 



rriy 

r]v + 2Lg—^ rjvpp 
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Appendix C : Loop integrals 

Here we collect the functions appearing in the loop contributions to the magnetic 
amplitude. 

The function K[q'^,mf,'m'j] appearing in the loop contributions to the magnetic am- 
plitude Mjy>zw Eq. (HD and M^^^ Eq. (|5|) is defined through 

/ 77rTDTr2 2TT77 ^ 2T = t9^''q'B,2[q',mlm^] + ... , (C.l) 

J {2tt) [r - m^J [[i - q) - m -\ 

as 2 

g^5,,[g^m^m?] |^ = (^) ir[g^m^m|] , (C.2) 



where in the MS scheme the term proportional to Aoo = 2) _ 4 ~ (^'(1) + ln(47r) + 1) 
appearing in B22[q'^ , m"^ , m^] has been subtracted. 

The F[q'^,m'^] function defined in the same loop amplitudes is 



with X = q^ jw? (note that this function is usually quoted in the literature with an 
extra factor of m?). 

For equal masses the K and F functions are related as 



6 - ) In + ^ 9^ + 4 m^] 



(C.4) 
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